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Abstract 

We study the symmetric 2-zero texture of the neutrino mass matrix, which is obtained 
from the symmetric Dirac neutrino mass matrix with 2-zeros and right-handed Majo- 
rana neutrino mass matrix with the general form via the seesaw mechanism, for the 
SUSY 50(10) GUT model including the Pati-Salam symmetry. We show that the 
only one texture in our model, having degenerate mass spectrum for the 1st and 2nd 
generation of right-handed Majorana neutrino, can simultaneously explain the current 
neutrino experimental data, lepton flavor violating processes and baryon asymmetry 
of the Universe. Within such a framework, the predicted values of the light and heavy 
Majorana neutrino masses, together with |?7e3|) <^CP ^-nd |(Mee)|, are almost uniquely 
determined. 
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1 Introduction 



We have now common information of neutrino masses and mixings [1]. Most remarkable one 
is that atmospheric neutrino mixing angle is almost maximal [2], while the solar neutrino 
mixing angle is large but not maximal [3, 4, 5], and further the ratio Aml^^/ Aml^^ is ~ 
with A ~ 0.2, which is much different from the quark mass spectra. Having estabhshed 
such precise measurements of dominant neutrino oscillation parameters, the maximal-large 
mixing angles with mass hierarchy of order A, we are now at a new stage of neutrino study. 
Our main concern is, not only how to reproduce the maximal-large mixing angles with less 
mass hierarchy but also how to predict the CP violating phases as well as the small f/gs [6]- 

Although there still remain many parameters of the neutrino mass matrix, we have al- 
ready grasped its global structure. Thus, it is an important task to exhaust the candidates 
for the models which are compatible with the present neutrino data, producing quite nat- 
urally the maximal-large mixing angles as well as the mild hierarchical mass ratio, and to 
make very strict predictions of all the neutrino parameters. Then, we can make the criterion 
of realistic models clear and such models may be checked and selected by the near-future 
experiments. 

It is well known that the following option for Mi with M^ (Georgi-Jarlskog type [7]) 



/ fld 

Ml = 



aa -3bd | m^. Ma = 

V 1 



/ arf " 

ad bd I rUb, (1.1) 
Vo 1, 



can reproduce the beautiful relations between the down-quarks and charged leptons at the 
GUT scale (Georgi-Jarlskog relations), 

TTir = rub, = 3ms, rUe = (1.2) 

with the (1-2) mixing as tan^^^^ = \/^ ^ ^ i^l- These relations are reahzed if we assume 
that each element of the down-quark Yukawa coupling is dominated by the contribution from 
either the 5 (10) or 45 (126) Higgs fields in the SU{5) (-SO(IO)) GUT, as follows: 

5(10) 

Ma, Mi; \ 5(10) 45(126) | . (1.3) 

Note that we need the "texture zeros", namely, some entries (the 1-1, 1-3 and 2-3 entries 
in this case) are far smaller than what we expect from naive hierarchical order of magni- 
tudes. Many people have studied the "texture zeros" extensively for the quark masses and 
mixings [8], and recently for the neutrino masses and mixings [9, 10]. It has been also con- 
sidered under the framework of the various GUTs [7, 11, 12, 13, 14, 15, 16], as we can see 
above. Such zero textures are most popular and may be some indication of family symmetry. 

Encouraged by the above fact, we further try to examine whether such simple assumption 
can work in the up-quark and neutrino sectors. In this paper, we adopt the so-caUed sym- 
metric four-zero texture [9, 17, 18] within the SUSY SO{10) GUT including the Pati-Salam 
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symmetry, which can relate not only to Mi, but also to M,^^. As for and Mi, in 
the symmetric four-zero texture, the following Higgs configuration [14] ^ 



can also reahze the Georgi-Jarlskog relation in eq. (1.2). In applying our GUT model, we 
need the information of the neutrino mass matrix with two zeros M,^ which is not directly 
obtained from M„, since it is related with the Dirac neutrino mass matrix with two zeros 
M^^ via the seesaw mechanism, M^ = Mj^M^^M^^, and we therefore have some freedom 
coming from the right-handed neutrino mass matrix Mji, to which only the 126 Higgs field 
couples, in order to determine which configuration of the Higgs representations for M^, M^,^ 
should be chosen to give proper neutrino masses and mixings. 

On the other hand, if the Nature demands the heavy Majorana masses for right-handed 
neutrinos to explain naturally the tiny neutrino masses via seesaw mechanism, the baryon 
number in the Universe may be affected by the leptogenesis which is caused by such heavy 
right-handed neutrino decay. Indeed the right-handed neutrino mass matrix plays a very 
important role in leptogenesis and it is considered one of the most hopeful scenarios to 
explain the origin of baryon number in the Universe, where the CP phases of the right- 
handed sector of neutrinos is very important. Combining the above information, what we 
should do next is to make definite predictions of various types of models which are compatible 
with the present experimental data, and to see what would be expected by including CP 
phases. In order to perform this, it is not enough to discuss the order of magnitude and 
we should make precise predictions based on strict theoretical arguments. Therefore, our 
senario may be one of the most hopeful approaches to make comparison of their predictions 
as definite as possible. 

In the previous papers [19, 20], we have shown that the symmetric two- zero texture of 
quark mass matrices can reproduce the neutrino maximal-large mixing angles by connecting 
them to lepton mass matrices by the Pati-Salam symmetry, with the right-handed Majorana 
mass matrix with four zeros. There the group coefficient factors are important to reproduce 
current neutrino experimental data. In this paper, we make a full analysis of such scenario in 
the SUSY SO{lQ) CUT and see how they are consistent with the neutrino masses and mixing 
angles as well as the baryon number in the Universe via leptogenesis, where the simplest 
form of right-handed neutrino mass matrix is extended to more general cases within 2- zero 
texture. Note an interesting fact that the original simplest form predicts two hghtest right- 
handed Majorana neutrino masses are degenerate. This is quite preferable if we want to 
explain the barion number generation of the Universe from the leptogenesis. At present we 
do not address what is the origin of these zero texture, leaving such more intersting question 
to the future task, which may be beyond the scope of this paper. 

This paper is organized as follows. In the section 2, the numerical analyses of masses and 
mixings are presented in the symmetric neutrino mass matrix with two zeros for the possible 
four textures of Mr. In sections 3 and 4, the lepton flavor violations and the leptogenesis 
are discussed in our model. Section 5 is devoted to summary. 

^In Ref. [14], Achiman and Greiner used this configuration in the five-zero texture. Our model is different 
from their model in this point. 



/ 10 
Md,Mi; 10 126 10 

V 10 10 




(1.4) 
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2 Symmetric two-zero texture in neutrino mass matrix 



2.1 The simplest form for Mr 



First, let us consider the following model-independent symmetric 2-zero texture including 
the CP violating phases, which we have investigated previously [20]; 





/O 




0^ 




/O 










a 
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rriy = 






I) 






h 








h 





(2.1) 



where j3 ~ 0(A), a ~ 0(1), — 0{1) and complex numbers, a^j3 and /i are converted to 
positive real numbers, a, /3, /i, by factoring out the phases with the diagonal phase matrix 
^. In such a texture, we examine how the parameters appearing in eq. (2.1) at the GUT scale 
are generally constrained from the present neutrino experimental data of sin^ 2^atm, tan^ ^sun 
and the ratio of Amg^^ to Am^t^^. As usual, we define the neutrino mixing angles which are 
expressed in terms of the MNS matrix [21]; 

Vmns = U}U,, 
where Ui and Uy diagonalizes Mi and M^,, respectively. 



diag(me, m^,m^), 
diag(m^,,m^^,m^J. 



(2.2) 



(2.3) 
(2.4) 



To examine the MNS matrix, we must take account of the contributions from the charged 
lepton side, Ui in eq. (2.2). The complex symmetric charged lepton mass matrix with 2-zeros 
is assumed to be written in terms of the real symmetric matrix M; ^ 

Ml = PiMiPi, (2.6) 

where Mi is diagonalized to Twf by real orthogonal matrix Oi [22] , 

OfMiOi ^ Mf'^-, (2.7) 

^This kind of 4-zero case has been studied extensively for the quark masses; 

/Oa„0\ /OadO\ 
Mu = flu by, Cu \ mt , Md= \ ad bd Cd \ nib . 
V c„ 1 / V Cd 1 / 

In this paper, the quark and lepton mass matrices are assumed to be factored out all the phases by the 
diagonal phase matrices in the 4-zero texture case. This is exactly possible in the case of 6-zero texture. 
Note that, however, we cannot factor out all the phases to make the matrix elements of M all real and there 
remains one phase as is seen in eq. (2.1). See Appendix A. 

^Here, we take the following symmetric matrix with 2-zeros for Mi, 

(0 ^in.emi_, \ 
^mefrifj, nifj, y^m^rn^ . (2.5) 
^JrrQm^ vrir j 

In general, one CP phase remains in this mass matrix, but its effect can be neglected as shown in Appendix 
B. 
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and, therefore, Mi is diagonalized by Pi*Oi as follows: 



■diag. 



(2.8) 



Similarly, it is supposed that the Dirac and right-handed Majorana neutrino mass matrices 
with two and four zeros are factored out the phases with the diagonal phase matrix P,^^ and 
Pr, respectively: 



M, 



R 



PrMrPr. 



(2.9) 
(2.10) 



On the basis where the charged lepton mass matrix is diagonalized, the neutrino mass matrix 
at Mr scale is obtained from cq. (2.1) 



M,{Mr) = Of Q' M,(Mr)QOi , 



(2.11) 



where 



M,(M, 



Q = PuPi* 




(2.12) 



(2.13) 



In order to compare our calculations with experimental results, we need the neutrino mass 
matrix at Mz scale, which is obtained from the following one-loop RGE's relation between 
the neutrino mass matrices at Mz and Mr [23] ; 



M,{Mz 









1 



1-tu 



M,{Mr) 







0" 




l-ea 



(2.14) 



















Here M^, is the neutrino mass matrix on the basis where charged lepton matrix is diagonalized 
(see eq. (2.11)). The renormalization factors and depend on the ratio of VEV's, tan/3^. 
Here, we ignore the RGE effect from Mqut to Mr scale considering that it almost does 
not change the values of masses for quarks and leptons. Using the form of eq. (2.14), we 
search the region of the parameter set {a, /3, h, (p, a, p) which are allowed by experimental 
data within 3(7 [1] : 

0.82 < sin^ 2e^tm , 
0.28 < tan^ ^sun < 0.64 , 



0.73 X 10 



< Am^,^ < 3.8 X 10-W 



— - cLim — 

5.4 X 10-^ < Am^„„ < 9.5 x lO^^eV^ . 



(2.15) 



In our model, wc show that these two large mixing angles can be derived from the 
symmetric four-zero texture with the Pati-Salam symmetry. We assume the following tex- 
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tures for up- and down-type quark mass matrices at the GUT scale [18], 



mb 
ma 







m.s 



mi-rrid 





nib -ma 

rrib - rrid 





mb 





rub 
mb_ 
V 



\ 



Irud 
mb 

1 



(2.16) 



which reproduces beautifully the down quark masses. For the up quark mass matrix, which 
is related to the Dirac neutrino mass matrix, we also take the following form 







mt 





mt 
rUc 

mt 
frnZ 
mt 



\ 



mi, 
mt 

1 



mt 



(2.17) 



This, together with the form of eq. (2.16), reproduces all the observed quark masses as 
well as CKM mixing angles. As for Md, it is well known that each element of M„ and is 
dominated by the contribution either from 10 or 126 Higgs fields, where the ratio of Yukawa 
couplings of charged lepton to down quark are 1 or —3, respectively. More concretely, the 
following option for (Georgi-Jarlskog type [14]) 



/ 10 
10 126 

V 10 



is known to reproduce very beautifully all the experimental data of mr^m 



rrif 



(2.18) 



as well as 



mi,,ras,md. On the other hand, M^, Mj,^ is related to M^^^, which is not directly connected 
to neutrino experiments, and we have not yet determined which configuration of the Higgs 
representations should be chosen to give proper neutrino masses and mixings. There are 16 
types of textures for M^, which are listed in Table 1. Once we fix their types, the Dirac 
neutrino mass matrix is automatically determined as 



M, 



/ 



^mumc 



mt 
. mc 



\ 



V 



mt 





rtit 



J 






(0 


a 




m = 
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b 






[0 


c 





with the Clebsch-Gordan (CG) coefficients denoted by *, which are either 1 or —3, according 
to whether the Higgs representation is 10 or 126. For the right-handed Majorana mass 
matrix, to which only the 126 Higgs field couples, we assume the following simplest texture: 





/o 


r 







r 










lo 





1 



(2.20) 
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Type 



Texture 



Type 



Texture 



126 
Si ( 126 10 10 
10 126 



^2 






126 





126 


10 


10 





10 


10 





126 






126 
Ai I 126 126 126 
126 126 



126 126 126 
126 10 



10 
10 10 10 
10 126 



A, 



10 
10 10 10 
10 10 



10 
Bi I 10 126 126 
126 126 



Bo 



10 
10 126 126 
126 10 



126 
Ci ( 126 10 126 
126 126 




126 
10 126 






126 


10 





10 





10 


10 


126 





126 


10 





126 






C2 



10 
10 10 126 
126 126 



C3 



126 
Fi ( 126 126 10 
10 126 



126 126 10 
10 10 



F2 



10 
10 126 10 
10 126 



F3 



10 
10 126 10 
10 10 



Table 1: Classification of the up-type mass matrices, and Mj 
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with two real parameters M3 and r. The neutrino mass matrix is now straightforwardly 
calculated as 



/ ^ 



2 



VO c(f + l) 




where a, b, c and d are defined in eq. (2.19). We know that the orders of the parameters in 
eq. (2.21) satisfy a <^ 6 ~ c <^ 1. In order to get a large mixing angle 623, the first term of 
the 2-3 element of M,y in eq. (2.21) should be of order d^, namely ac/r ~ 0{d'^). This fixes 
the value of r as 

-10-(^-^\ (2.22) 

which is indeed the ratio of the the right-handed Majorana mass of the 3rd generation, M3, 
to those of the 1st and 2nd generations, Mi and M2. With this smaU r, M,y is approximately 
given by 

/O — 0\ 2 /'J/0'U\,2 2 



with 






r 







2ab 


ac 
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r 





ac 
r 


rf2 




ac 2a6 ^ a 



2 



where (3 <^ a and /i ~ C'(l). The forms of /i, a and /3 are written in terms of rrit, rric and 
m„, with a parameter r, or equivalently h. In Table 2, wc can classify the 16 types into five 
classes S, A, B, C and F . We have shown that the types belonging to a corresponding class 
yield the same predictions for mixing angles and masses [19]. 

Now, we can predict the values of a and j3 from the up-quark masses at the GUT scale, 

ruu = 0.36 ~ 1.28 MeV, (2.25) 
ruc = 209 ~ 300 MeV, (2.26) 
mt = 88 ~ 118 GeV, (2.27) 

which are obtained taking account of RGE's efi'ect to the quark masses at the EW scale [8] . 
We have shown that the allowed region of a and /? given by a neutrino mass matrix with 
two zeros in eq. (2.1) and the region of a and (3 for the five classes in our model, which 
are predicted from the up-quark masses at the GUT scale, are slightly separated on the 
a-(3 plane, as seen in Figure 1. However, the light quark masses are ambiguous because 
of the non-perturbative QCD effect. Therefore, the allowed region of up quark mass, m„, 
may be enlarged ^ . We obtained the overlapped region for the type Si (of the class S), 
as seen in Figure 4 of Ref. [20], enlarging the values of up quark mass at the GUT scale, 
m„ = 0.36 ~ 2.56 MeV. In the overlapped region, the allowed values of the parameters 
including the phases are restricted to very narrow regions. On the other hand, our results are 

®Thc other possibility, that is, the effect of deviation from rric in the 2-2 element of M„, are discussed in 
Appendix C. 
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Type 


(f 


r in ac/h unit 


a in 2mc/ yJmuTnt unit 


[3 m -^mc/mt unit 




9 


-1/3 


1 


-3 




1 


—3 


1 


—3 


^9 


1 


1 


1 


1 


A2 


1 


9 


1 


1 




9 


1/9 


1 


1 


A I 


1 


1 


1 


1 


Bi 


9 


-1/3 


1 


-1/3 


B2 


1 


-3 


1 


-1/3 


Ci 


9 


1 


-1/3 


1 


C'2 


!) 


-1/3 


-1/3 


-1/3 


C3 


1 


-3 


-1/3 


-1/3 


C4 


1 


9 


-1/3 


1/9 


Fi 


9 


-1/3 


-3 


-3 


F2 


9 


1/9 


-3 


9 


F3 


1 


1 


-3 


1 


Fi 


1 


-3 


-3 


-3 



Table 2: The forms of h, a and (3 with the value of d for each type. 



Type 


sin^ 023 


tan^ 012 


Ami2 (eV^) 


Amil (eV^) 


m^i (eV) 


(eV) 


m^j (eV) 




Si 


0.42 ~ 0.43 


0.28 


3.8 X 10-3 


5.4 X 10-5 


0.0014 


0.0075 


0.062 


TT TT 


S2 


0.43 


0.28 


3.8 X 10-3 


5.4 X 10-5 


0.0014 


0.0075 


0.062 





Type 




Jcp 


l(Mee)| 


Ml (GeV) 


M2 (GeV) 


M3 (GeV) 


P 


Si 


0.010 ~ 0.025 


(-4.8 ~ -1.1) X 10-3 


0.019 ~ 0.024 


1.1 X 10'' 


-1.1 X 10*^ 


3.0 X 10^^ 


- ~7r 


S2 


0.010 





0.024 


1.1 X 10*^ 


-1.1 X 10" 


3.5 X 10" 


TT 



Table 3: The allowed values for the neutrino observable in the types S'l and 5*2. Mi, M2 
and M3 are the masses for the 1st, 2nd and 3rd generation of the right-handed Major ana 
neutrino. 

almost independent of the phase parameter a and therefore we take cr = in our calculations 
for simplicity. By taking those values of parameters, we can obtain the prediction of \Ue3\, 
Jcp, \{fTiee) \ and the absolute masses of neutrinos, as shown in the Table 3. In the Table 3, 
we also list the allowed values for sin^^23) tan^^^i2, Am|2. ^77132, Mi, M2 and M3 in the 
types Si and S2. As we can see that, from the Table 3, the difference between the types S'l 
and 5*2 is only the scale of M3; the order of M3 for Si is larger 1 order than 5*2. For the type 
S'2, we obtained the overlapped region, enlarging = 0.36 ~ 2.64 MeV. Therefore, the 
type S'l is the best type in the class S. 

In the next section, we will examine whether or not more general form of the right-handed 

Majorana neutrino mass matrix, which leads to the neutrino mass matrix with two zeros 
of eq. (2.1), can have parameter regions consistent with the neutrino experimental data, 
without enlarging the values of up quark mass at the GUT scale. 
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2.2 Including new parameters in Mr 



2.2.1 Properties for the general form of Mr 

Up to here, we have assumed the simplest form for Mr, having two parameters, 



Case I : M, 



R 




(2.28) 



However, the actual case might have a general form which leads to the neutrino mass matrix 
with two zeros in eq. (2.1). Thus, we here include the new parameters s and t in Mr as 
follows: 



Mr^ 



/O r 0\ 

r s t 
\0 t l) 



Ms, 



(2.29) 



where r, s and t are taken to be real. In order to clarify each effect of the new parameters, 
let us examine the above form, eq. (2.29) by dividing it into the following three cases: 







(0 


r 






Case II : 


Mr^ 


r 


s 


I] 


Ms, 






lo 













fo 


r 






Case III : 


Mr^ 


r 





t 


Ms, 








t 












r 






Case IV : 


Mr = 


r 


s 




Ms. 






u 


t 







(2.30) 



(2.31) 



(2.32) 



The final expression of the left-handed Majorana neutrino mass matrices for each case can 
be obtained via seesaw mechanism: 



Case I : 



Case II : 



Case IV : 



Case III : M^ 



/ ^ 



a 
r 



V 
/ 



a 
r 



Vo 

/ 

of. 

r 

V 

/ 

r 

V 



2ab I ^2 



2ab 




Ms' 



a 

^- ^t + cd 

r r 

r 

2ab _ _ 2ac^ _^ a^^2 _^ ^2 





-^t + cd 
d^ 



m 



r 
ad 



t + cd 





^-^t + cd 

r r 
d^ 



Ms 



(2.33) 



(2.34) 



(2.36) 
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As is seen from the above expressions, the parameter s affects only on the 2-2 element of 
M^, and the additional contribution to the 2-3 element comes only from the parameter t. 
In eqs. (2.33) and (2.34), we found that the 2-3 element of M^y in the case 11 is the same as 
that in the case 1, namely, the condition for getting a large mixing angle 623 in the case 11 is 
the same as that in the case 1. On the other hand, the condition for the case 111 is, similarly, 
the same as that in the case IV. Note that the 1-2 element of Mi, has the same form in all 
case. 

For each case, we will search for the parameter regions, which are allowed by the experi- 
mental data within 3a, and show that, in the case II, the types Si and S2 have the allowed 
regions, the types Ci, F3 and F4 in the case III, and the types Si, S2, Ai, Bi, Ci, F3 and F4 
in the case IV. 



2.2.2 Case II 

First, let us examine the case II. Since a new parameter, s, affects only on the 2-2 element 

of Mj,, namely, a, the allowed regions for the classes S, A, B, C and F in the case I, as 
depicted in Figure 1, can be enlarged only in the direction of a. The first term in the 2-2 
element becomes comparable with the second term at the following value of s: 

2ab a?s 2br , q ^ , „ 7 „„x 
~ — — ^A'^-Sx 10-^. (2.37) 

Thus, the difference between the region of a in the case I and the case II is appreciable 
around this value of s. We obtained the overlapped region in the types Si and 5*2 with 
h— 1.3, as depicted in Figure 2 and 3. It is shown that the region which is consistent with 
the experimental data has now been focused only on the narrow region. The typical values 
for these types arc listed in Table 4, where the values at the maximum and minimum value 
of If/esI are written. In this table, as we have expected, we find that the overlapped regions 
can be obtained around s ~ 10~^^~^''. The predicted values of \Ue3\ is given 

I C/e3| = 0.014- 0.071, (2.38) 

for the type 5"!, and 

It/eal = 0.013-0.016, (2.39) 

for the type S2. 



2.2.3 Case III 

Next, we discuss the case III, in which the form of the 2-3 clement of M^y is different from the 
one in the case I. We have listed the condition for getting a large 623 of each class in Table 5. 
Here, a new parameter t is included in the condition of r. Therefore, the allowed region for 
five classes, as seen in Figure 1, can be enlarged both in direction of a and f3. The first term 
in the 2-3 element becomes comparable with the second term at the following value of t: 

— ~— t^t~-~A^~2x 10"^ 2.40 
r r a 
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Type 


^max 


^min 


^max 


^miii 




0.071 


0.014 


0.016 


0.013 


</> 


0.0 


0.0 


0.0 


0.0 


P 


57r/8 


TT 


TT 


TT 


l(Mee)| 


0.022 


0.027 


0.029 


0.027 




-0.01 


0.0 


0.0 


0.0 


m,.! (eV) 


1.7 X 10~^ 


1.4 X 10"^ 


1.4 X 10-3 


1.4 X 10-3 


rrijy^ (eV) 


7.9 X 10 


7.5 X lO""^ 


7.8 X 10 


7.5 X 10 


TO,^3 (CV) 


5.8 X 10 


5.8 X 10 


5.8 X 10 


5.9 X 10 


7—— 

sill" 6*2:! 


0.4.") 


0.43 


0.44 


0.43 


tan^ 6*12 


0.29 


0.29 


0.28 


0.28 


Ami2 (eV") 


3.3 X 10-3 


3.3 X 10-3 


3.3 X 10-3 


3.4 X 10-3 


Ami, (eV") 


6.0 X 10"^ 


5.4 X 10"^ 


5.9 X 10-^ 


5.5 X 10-5 


m„ (McV) 


0.56 


0.52 


1.04 


1.00 


mc (MeV) 


300 


230 


300 


280 


TOt (GeV) 


88 


88 


108 


108 


s 


1.0 X lO"*^ 


6.5 X 10-^ 


7.0 X 10-" 


6.0 X 10-'^ 


Ml (GcV) 


-2.7 X 10^ 


-2.7 X 10^ 


-9.6 X 10^ 


-9.6 X 10^ 


M2 (GeV) 


3.0 X lO'-* 


2.0 X 10^ 


3.6 X 10^ 


3.1 X 10^ 


M3 (GeV) 


3.0 X 10^^ 


3.0 X 10^^ 


5.0 X 10" 


5.0 X 10" 



Table 4: Results of numerical calculations for 5*1 and 5*2 in the case II. The values at the 
maximum and minimum value of |t/e3| are written. 



Thus, the difference between the region of a and /5 in the case I and the case III is appreciable 
around this value of t. We show that the types Ci, F3 and F4 provide the overlapped region 
as depicted in Figure 4, 5 and 6. The typical values of Ci, F3 and F4 are listed in Table 6. 
In similar to the case II, as we have expected, we find that the overlapped regions can be 
obtained around t ~ 10~^. Then, the predicted values of \Ue3\ is given for the type Ci, 

l^esl - 0.025, (2.41) 

for the type F3, 

l^/esl = 0.059- 0.17, (2.42) 



and for the type F4, 



\Ue3\ ^ 0.19. (2.43) 



The types Ci and F4 have very narrow region. On the other hand, wide region is obtained 
in the type F3. 



2.2.4 Case IV 

Finally, we consider the most general case for Mr, which includes two new parameters, s 
and t. In this case, because eq. (2.36) includes both cqs. (2.34) and (2.35), we can expect 
that 81,82, which are allowed in the case II, and Ci,F3,F4, which are allowed in the case 
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Type 


condition 


Type 


condition 


^1 


», -a(,c+3tj 
— 3/i+c 




h—c 




h—c 




^, 3a(3c+i) 
— /i+3c 


A 

A3 


~ 3(3fe+3c) 


/I 


a(c—t) 

rp 1-^ 1 \ 1 

h—c 


Bi 


^ -a(c-i) 
' — 3(/i-c) 


B2 


^ -a(3c+t) 
~ /i+3c 


Ci 


a(c— t) 
h—c 


C2 


„ ^. -a{c-t) 
' — 3{h-c) 


C3 


-a(3c+t) 
~ /i+3c 


C4 


3a(3c+t) 
~ h+3c 




-a(c+3t) 
~" 3/i+c 




~ 3(3/i+c) 




a(c— t) 

/p ^ 1 V y 
/l— C 


F4 


^ ^ , -3a{c-t) 
h—c 



Table 5: The condition for getting a large 623 in the case III and IV. 



Type 


Ci 




jpmni 


F4 




0.025 


0.17 


0.059 


0.19 





0.0 


-7r/8 


7r/8 




p 


7r/8 


7r/2 


7r/8 


37r/4 


mee)\ 


0.019 


0.027 


0.084 


0.57 


Jcp 


0.0042 


0.033 


-0.005 


0.031 


m^i (eV) 


1.4 X 10-^ 


3.1 X 10-3 


1.6 X 10-3 


5.7 X 10-3 


(eV) 


7.5 X 10"^ 


9.6 X 10-3 


8.3 X 10-3 


9.9 X 10-3 


m^3 (cV) 


6.2 X 10"^ 


4.4 X 10-^ 


3.7 X 10-^ 


4.4 X 10"^ 


sin"^ 6*23 


0.46 


0.45 


0.49 


0.30 


tan^ 012 


0.28 


0.40 


0.29 


0.56 


Amis (eV^) 


3.8 X 10-3 


1.9 X 10-3 


1.3 X 10-3 


1.8 X 10-3 


Amii (eV^) 


5.5 X 10"^ 


8.2 X 10-^ 


6.7 X 10"^ 


6.6 X 10-5 


TUu (MeV) 


0.48 


1.2 


1.24 


1.12 


nic (MeV) 


290 


300 


220 


220 


mt (GeV) 


92 


93 


88 


113 


t 


1.4 X 10-3 


3.0 X 10-3 


3.0 X 10-3 


2.5 X 10-3 


Ml (GeV) 


1.2 X 10' 


4.7 X 10'^ 


6.5 X lO'' 


6.2 X 10"^ 


M2 (GeV) 


-6.1 X lO'J 


-4.5 X 10'' 


-4.5 X 10" 


-5.6 X 10" 


Ms (GeV) 


3.09 X 10^^ 


5.0 X 10" 


5.0 X 10" 


9.0 X 10" 



Table 6: Results of numerical calculations for Ci,F3 and F4 in the case III. The values at 
the maximum and minimum value of |C/e3| ^-re written. 



12 



Type 


^max 


^min 


^max 


^min 




0.027 


0.014 


0.19 


0.021 


<i> 


0.0 


0.0 


-7r/8 


0.0 


P 


77r/8 


TT 


7r/4 


0.0 


\{Mee)\ 


0.024 


0.027 


0.27 


0.034 




-4.0 X 10"^ 


0.0 


0.034 


0.0 


m^i (eV) 


1.5 X 10"^ 


1.4 X 10"'^ 


4.3 X 10-3 


1.7 X 10-3 




7.6 X 10"'^ 


7.6 X lO""^ 


9.1 X 10 


7.5 X 10-3 


TO,^3 (CV) 


5.8 X 10 


5.8 X 10 


5.0 X 10 


5.9 X 10 


sill" 02:i 


0.44 


0.43 


0.38 


0.4(i 


tan^ 012 


0.29 


0.28 


0.39 


0.35 


Ami2 (eV") 


3.3 X 10-3 


3.3 X 10-3 


2.4 X 10-3 


3.4 X 10-3 


Ami, (eV") 


5.5 X 10"^ 


5.5 X 10"^ 


6.4 X 10-5 


5.4 X 10-5 


m„ (McV) 


1.28 


0.52 


0.96 


0.48 


TOc (MeV) 


280 


290 


270 


280 


mt (GeV) 


88 


88 


93 


108 


s 


1.0 X 10-*^ 


1.0 X 10-" 


8.0 X 10"'^ 


6.0 X 10-** 


t 


1.0 X 10-"* 


1.0 X lO-'^ 


2.0 X 10^3 


4.0 X 10-3 


Ml (GeV) 


-1.5 X 10» 


-2.9 X 10^ 


-2.9 X 10^ 


1.1 X 10^ 


M2 (GeV) 


3.1 X 10" 


3.0 X 10" 


2.0 X 10" 


-5.0 X 10" 


Ms (GeV) 


3.0 X lO'*' 


3.0 X lO'*' 


5.0 X 10'* 


5.0 X 10'* 



Table 7: Results of numerical calculations for and 5*2 in the case IV. The values at the 
maximum and minimum value of |C/e3| ^-re written. 



Ill, have the overlapped region. By numerical calculations, we have confirmed that Ai and 
Bi also provide the overlapped regions as depicted in Figures 7 ~ 13. The typical values of 
these seven types are listed in Table 7, 8 and 9. We present the predicted values of {Uesl for 
seven types as follows: 



\Ue3\ 


= 0.014- 0.027 (Type Si), 


(2.44) 


lUesl 


= 0.021 - 0.19 (Type ^2), 


(2.45) 


\Ue3\ 


= 0.032- 0.17 (Type Ai), 


(2.46) 


\Ue3\ 


= 0.049 - 0.11 (Type 


(2.47) 


\Ue3\ 


~ 0.012 (Type Ci), 


(2.48) 


\Ue3\ 


= 0.056- 0.19 (Type F3), 


(2.49) 


\Ue3\ 


~ 0.17 (TypeF4). 


(2.50) 



In conclusion, the prediction of |?7e3| depends on the types of Dirac neutrino mass matrix 
and right-handed Majorana mass matrix considerably. 

3 Lepton flavor violations : ej 6^7 

In the model of MSSM with right-handed neutrinos, lepton flavor violations (LFV) are 
induced through the renormalization group effects to the slepton mixings and the predicted 
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Type 


^max 




^max 


^min 




0.17 


0.032 


0.11 


0.049 


</> 


0.0 


0.0 


0.0 


0.0 


P 


TT 


0.0 


tt/2 


0.0 


l(Mee)| 


0.24 


0.049 


0.096 


0.075 


-^CP 


-0.01 


0.0 


0.012 


0.0 


m,.! (eV) 


3.6 X 10"^ 


1.6 X 10-3 


1.9 X 10-3 


1.9 X 10-3 


rrijy^ (cV) 


8.6 X 10 


8.4 X 10 


8.5 X 10 


8.5 X 10 


TO,^3 (cV) 


6.0 X 10 


5.1 X 10-^ 


4.7 X 10 


4.7 X 10 


■) 

sill" O-z.i 


0.42 


0.48 


0.48 


0.4!) 


tan^ 6*12 


0.40 


0.31 


0.29 


0.37 


Ami2 (eV") 


3.5 X 10-3 


2.6 X 10-3 


2.1 X 10-3 


2.1 X 10-3 


Ami, (eV") 


6.2 X 10"^ 


6.7 X 10-^ 


6.8 X 10-^ 


6.8 X 10-5 


m„ (McV) 


0.36 


0.44 


0.8 


0.8 


mc (MeV) 


210 


260 


300 


300 


rrit (GeV) 


113 


108 


103 


103 


s 


3.0 X 10-** 


4.0 X 10-" 


8.0 X 10-" 


8.0 X 10-" 




1.5 X 10-3 


1.5 X 10-3 


3.0 X 10-3 


3.0 X 10-3 


Ml (GeV) 


-1.9 X 10'^ 


-4.4 X lO'' 


3.5 X 10" 


3.5 X 10" 


M2 (GeV) 


3.8 X 10" 


8.8 X 10*^ 


-5.0 X 10*^ 


-5.0 X 10*^ 


Ms (GeV) 


5.0 X lO'*' 


5.0 X lO'*' 


5.0 X 10'" 


5.0 X 10'" 



Table 8: Results of numerical calculations for Ai and Bi in the case IV. The values at the 
maximum and minimum value of |C/e3| are written. 



branching ratios of the processes can be comparable with the current experimental upper 
bound [24, 25, 26, 27, 28, 29, 30] : 

Br(^ ^ 67) < 1.2 X 10"" [31] , (3.1) 
Br(T ^ 67) < 3.6 X 10-^ [32] , (3.2) 
Br(T ^ /X7) < 3.1 X 10"^ [33] . (3.3) 

Therefore, we have to examine these decay rates carefully in our model. 

Let us start with writing down the leptonic parts of soft SUSY breaking terms as follows: 
-£soft = {mDijLlLj + {ml)ije*j^^Rj + {ml)iji>*i^i>Rj + {Ae)ijHde*j^Lj + {A^)ijH^D*i^Lj ,(3.4) 

where Lj, CRi and VRi are the supersymmetric scalar partner of left-handed lepton doublet, 
right-handed lepton singlet and right-handed neutrino, {m'j^)ij, {ml)ij and {ml)ij are the 3x3 
hermitian slepton mass matrices, {Ae)ij and {A^)ij are trilinear couplings of Higgs doublets 
Hu,Hci and sleptons, respectively (A-term). Hu and Hd couple to neutrinos and charged 
leptons, respectively. Non-vanishing off-diagonal elements of slepton mass matrices become 
the new source of LFV. 

In minimal-supergravity (mSUGRA) models, it is assumed that the slepton mass matrices 
are diagonal and have common mass scale mg at the GUT scale, and that the trilinear 
couplings are proportional to Yukawa couplings: 

/ 2\ / 2\ / 2\ ST 2 2 2 2 

At, = Y^aQUiQ , Ae = Yeaouio , (3.5) 
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Type 


Ci 


77' max 
^3 


TTimm 
^3 






\Uez\ 


0.012 


0.19 


0.056 


0.19 


0.17 


<i> 


0.0 


-7r/8 


0.0 


0.0 


-7r/8 


P 


0.0 


37r/4 


0.0 


TT 


77r/8 




0.025 


0.35 


0.089 


0.67 


0.51 




0.0 


0.034 


0.0 


0.0 


0.031 


m^i (eV) 


1.4 X 10~^ 


2.9 X 10-3 


1.6 X 10-3 


5.5 X 10^3 


5.6 X 10^3 


TTiy^ (eV) 


7.5 X 10""^ 


8.6 X 10 


7.5 X 10 ^ 


1.1 X 10-^ 


9.9 X 10 


m,^3 (cV) 


6.1 X 10^"^ 


4.0 X 10 


3.7 X 10 


3.9 X 10 


4.8 X 10 


sill" 62;', 


0.4() 


0.44 


O.-IO 


0.39 


0.30 


tan^ 012 


0.28 


0.35 


0.35 


0.56 


0.59 


Ami2 (eV^) 


3.7 X 10-3 


1.5 X 10-3 


1.3 X 10-3 


1.4 X 10-3 


2.2 X 10-3 


Amii (eV^) 


5.5 X 10-5 


6.6 X 10-5 


5.4 X 10-5 


9.4 X 10-5 


6.7 X 10-5 


m„ (MoV) 


0.88 


1.12 


1.24 


1.12 


1.2 


me (McV) 


300 


290 


210 


210 


210 


mt (GeV) 


88 


88 


88 


113 


118 


s 


3.2 X 10- ^ 


2.0 X 10- ^ 


4.0 X 10-'^ 


4.0 X 10- ' 


1.0 X 10- ' 


t 


1.9 X 10-3 


3.0 X 10-3 


3.0 X 10-3 


2.5 X 10-3 


2.5 X 10-3 


Ml (GeV) 


2.8 X 10^ 


4.5 X 10" 


6.5 X 10" 


6.3 X 10" 


6.7 X 10" 


M2 (GeV) 


-9.4 X 10*^ 


-4.4 X 10*^ 


-4.3 X 10*^ 


-5.3 X 10^ 


-5.5 X lO'* 


Ms (GeV) 


2.85 X 10'" 


5.0 X 10'* 


5.0 X 10'* 


9.0 X 10'* 


9.0 X 10'* 



Table 9: Results of numerical calculations for Ci,F3 and F4 in the case IV. The values at 
the maximum and minimum value of \Uez\ are written. 



and the same conditions are assumed in quark sector. The mass of supersymmetric fermion 
partner of gauge bosons (gauginos) are also fixed to be M1/2 at the GUT scale. Even if no 
sources of LFV are assumed at the GUT scale, the LFV will be induced in slepton mass 
matrix through renormalization of Yukawa and gauge interactions. The one-loop renormal- 
ization group equation (RGE) for left-handed slepton mass matrix is given by 



d , o . d 
1 



MSSM 



+ 



:m|yJ-y, + yJ-y,m|),,- + 2(yJ-m,y, + m^^yJy, + 4^,),,- , (3.6) 



where the first term is the MSSM term which is lepton flavor conserving, while the second 
term contains the source of LFV, Y,, is the neutrino Yukawa coupling matrix (= My^/vu). 

It is easy to see in eq. (3.6) that the Yukawa coupling of the neutrino contributes to the 
LFV. Assuming the boundary conditions of eq. (3.5), we obtain the leading log approximation 
for the off-diagonal elements of left-handed slepton mass matrix at the scale of right-handed 
neutrino masses as follows [25, 26, 30]: 

1 

Stt' 

where the matrix H^j is defined by 

H, . [Y^^LY:), , L ^ diag (in In In ^) , (3.8) 



K)«^-^-^(3^o + «o)^^. , (3-7) 



15 



where is the neutrino Yukawa couphng matrix on the basis where the charged lepton and 
right-handed Majorana mass matrix is diagonal. The branching ratio for the LFV processes: 



ej7 is approximately given by 



BrCe, 



6^7) - 7^ 



tan^/3„ , 



(3.9) 



where ms is the typical mass scale of superparticles, a ~ 1/128 and Gp is the Fermi coupling 
constant, respectively. The excellent approximation of ms to the exact RGE result is given 
by [34] 

m| ~ 0.5 mlM^/^{ml + 0.6M^/2f ■ (3-10) 

It is clear that the element H12, i?i3 and H23 dominantly contribute to the processes of 
// — > 67, T — > 67 and r — > /i7, respectively. 

Let us calculate the matrix Hij in our model. The is given by 

Z = 0nY,0i, (3.11) 

where matrix Oji and Oi are the orthogonal matrices which diagonalize the right-handed 
Majorana mass matrix and charged lepton mass matrix, respectively. Then, the matrix Hij 
is given by 



Hi, = {OjY}Ol)ikLk{ORY,Oi)kj . 



For the case I, matrix Or is given by 

Or 



( 1/v^ \ 

-1/V2 1/V2 

1/ 



(3.12) 



(3.13) 



then we obtain the formulae of Hij assuming the type Si for neutrino Yukawa coupling 
matrix which is the mostly allowed one : 



H12 ~ 



(OmO/12 + O121O122) + (bOi2i + cOisi) (bOi22 + cOm) 

Mgut 



a {h{0m0i2i + O111O122) + c{0ii20m + Cm 0^32)}] In ■ 



Ml 



M. 



GUT 



H 
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+(cOu2 + dOi^i){cOi22 + dOi^2) In ■ 

{OinOas + OmOm) + (6O21 + cO^ai) {hOm + cOm) 



M3 



(3.14) 



a {h{On^Oi2i + OniOm) + c(Ou30i3i + 0^0^)}] In 



M, 



GUT 



Ml 



M, 



GUT 



/^23 - 



+(cOi2i + dOi3i){cOi23 + dOm) In ■ 
a' (Ou20u3 + 0,220,23) + {hOi22 + cOz32) [hOm + cO^aa) 



M3 



(3.15) 



a {&(Ou30,22 + O112O123) + c{On30i32 + 0,120,33)}] In 



M 



GUT 



Ml 



+ (cO,22 + dOi32){cOi23 + (^0,33) In 



M, 



GUT 



M3 



(3.16) 



16 



where we have used the mass spectrum for the case I: Mi — M2. In the following calculations 
of branching ratios, the relations : a = — Sou, b = bu,c = Cu,d = —3du for the type 5*1 and 
a = — 3a„, h = hu,c = Cu,d = for the type 5*2 will be taken. The orthogonal matrix Oi is 
approximately given by [18] 



Oi 



( 



\ 



/ me 



\ 



I 




e = 



rrif, rrip , 



Under the above parameterization, we can calculate the off-diagonal elements of matrix if, 



tl2 



!13| 



'23| 





, , GUT 

ab m — — — 


-6^£ln 


Mgut 




Ml 


Ml 




, ^GUT 

acm — — — 


— hc£ In 


Mgut 




Ml 


Ml 




, 1 -/W^GUT 

cam — — — 


-c'51n 


Mgut 




M3 


Ml 



-dVsln 



M, 



GUT 



+ d^8e In 



M3 

Mgut 



c^d^^^sln^, (3.18) 
M3 ' ^ ' 



M, 



~ (i In 



M, 



gut 



Ms 



Mr 



M, 



Mr 



M, 



(3.19) 
(3.20) 



where the terms including the charged lepton mixing matrix are dominant ones. These 
formulae provide the following relations of branching ratios : 



Br(// 67) < Br(T 67) < Br(T n'j) 



(3.21) 



for the type Si and 32- For the case 11, 111 and IV, the predicted branching ratios are almost 
the same as the case 1. Numerical calculations of the branching ratios using the leading log 
approximation (eps. (3.7), (3.8), (3.9) and (3.10)) of — >• 67, r — >• 67 and r — > /i.7 are 
presented in Figure 14, 15 and 16, respectively ^. In these figures, the branching ratios 
are scatter plotted in the region of mo = 100 ~ 1000 GeV, M1/2 = 100 ~ 1000 GeV and 
tan /? = 5 ~ 50 for each processes. The Qq is fixed to be Qq — 0. 

As seen in these figures, the branching ratios of all processes are safely predicted below 
the current experimental upper bounds. The predicted branching ratios for the type 5*2 
are lower almost one order of magnitude than the one for the type 5*1. This is due to the 
differences in CG coefficients in 3-3 elements of neutrino Yukawa coupling matrix. Only 
T — s> /i7 process for the type ^i may be observed in the future experiments in which the 
sensitivity will reach to be Br(r ^^7) < 10~^ [35]. The predicted branching ratios of 
T ^ H'-f process for the type 5*2 and the other processes for both types are too small to be 
observed even in the future experiments. 



^Note that this approximation deviate significantly from exact RGE result in the region of large M1/2 
and small mo [34]. However, this deviations are at most of a factor ~ 10. For mo = 100 Gev, a discrepancy 
between full RG results and leading log approximation is of about one order of magnitude at M1/2 ^ 1 TeV, 
while for mo = 300 GeV, this is reduced to be about a factor of two. The size of discrepancy depends weakly 
on the scale of right-handed Majorana neutrino masses [34]. 
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4 Thermal Leptogenesis 



In this section, we discuss the calculation of baryon asymmetry of the universe based on the 
leptogenesis scenario [36] for our textures. In the leptogenesis scenario, lepton asymmetry 
is generated by the CP violating out-of-equilibrium decay of heavy right-handed Majorana 
neutrinos. Let us consider the CP asymmetry parameter ei, which is generated in the decay 
of i-th generation of right-handed Majorana neutrino Ni. The is defined as 

^ ^ TiN,^HL)~TiN,^HT) 
' r{Ni^ HL) + r{Ni^ HL) ' 

where H and L are the ordinary Higgs and lepton doublet. At tree level, the decay width of 
Ni can be easily calculated as : 

po ^ (nyj)n ^. _ (4 2) 

As seen in (4.2), even if the Yukawa coupling matrix Y^, contains complex elements, CP 
symmetry is not violated at tree level. Therefore, we should consider the one-loop contribu- 
tions. It is well-known that CP is violated in the interference between the tree diagram and 
one-loop self-energy and vertex correction diagrams. Summing up the one-loop vertex and 
self-energy corrections, the CP asymmetry is given by 

e^ = er + ef = -^—^EMiYMU [v (x) + s (x)] , (4.3) 

where v{x) and s{x) are the self-energy and vertex correction functions with x = M'j /Mf 
[37]. In the minimal supersymmetric standard model (MSSM) with right-handed neutrinos, 
they are given by [37, 38] 

vix) - J-x c(x)- (^^^1-^^)^-^^- - (^-l)v^ a 4) 

v[x) - m ^ J , s[x) _ _ Ml)2 + M2(r0)2 " (^ _ 1)2 + (r^.y/Mi ' ^^-^^ 

which is available for both cases of the hierarchical case and the quasi-degenerate case of Mi 
and Mj. 

In order to calculate the baryon asymmetry, we need to solve the Boltzmann equations in 
thermal leptogenesis scenario [39] . We can use the approximate solution of these Boltzmann 
equations as 

r]BC^0Mj2eiKi , (4.5) 

i 

where rjB is baryon asymmetry of the universe, Ki is so-called dilution factor which describe 
the wash-out effect of generated lepton asymmetry. The Ki is approximated as [40] 

f lO;hV \ ( \-'' ^ {M^DMlj,)ii 
Ki ~ 0.3 — In , mi = , ^ . (4.6) 

In the following, we compare the current range of observed baryon asymmetry [41] : 

77b = (6.2 - 6.9) X 10"^° . (4.7) 
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with the predicted values of our model. 

It is convenient to discuss the hierarchical case : Mi <^ M2, and degenerate case : 
Ml ~ M2 <S M3 of right-handed Majorana masses, separately. For the hierarchical case, 
from the model independent analyses of thermal leptogenesis [42, 43], the lightest Majorana 
neutrino mass must satisfy the condition : 



to generate the observed baryon asymmetry of the universe. The case II, III and IV corre- 
spond to the hierarchical case. In these shown in the numerical results of tables, 
the lightest Majorana neutrino mass Mi is lighter than 4.9 x 10^ GeV for all types to satisfy 
the conditions of the current neutrino experiments. Therefore, it is impossible to explain the 
observed baryon asymmetry by thermal leptogenesis for the case of II, III and IV. 

On the other hand, the case I corresponds to the degenerate case : Mi ~ M2 -C M3. 
This case satisfies the relation : M2 — Mi < Fi -|- F2. It is easy to find in eqs. (4.4) that there 
occurs an enhancement of CP asymmetry for some region of the degeneracy and = for 
the case of exact degeneracy : Mi = M2. The scenario utilizing this enhancement is called as 
"resonant leptogenesis" [44, 45]. Some author showed that observed baryon asymmetry can 
be generated with considerably light right-handed neutrino masses, in complete accordance 
with the current solar and atmospheric neutrino experiments [45, 46, 47, 48, 49, 50, 51, 52]. 
This is a candidate to solve the gravitino problem [53] . 

The mass eigenvalues of right-handed Majorana neutrinos for the case I are exactly 
degenerate : Mi = M2 = rMs ~ 10^ GeV. However, it is natural to explain that the 
mass spectrum may be somewhat deviated from exact degeneracy by, for example, quantum 
corrections. Let us define the degree of degeneracy for Mi and M2 by 



The predicted baryon asymmetry is shown in Figure 17 as a function of AM. As seen in 
Figure 17, if the degree of degeneracy is the level of AM ~ 10~^, the predicted asymmetry 
is consistent with the observed value for the type Si in the case I. It is concluded that the 
baryon asymmetry can be explained by the resonant leptogenesis scenario in the suitable 
region of the mass degeneracy in our model. Almost the same result is obtained for the type 



We have investigated the symmetric 2- zero texture of neutrino mass matrix for the possible 
four textures of the right-handed Majorana neutrino together with the Dirac neutrino mass 
matrix with two zeros, under the SUSY 5*0(10) GUT model including the Pati-Salam sym- 
metry. We made a full analysis for the parameters included in such four cases of neutrino 
mass matrices and showed how they arc consistently explain the neutrino masses and mixing 
angles as well as the baryon number in the Universe via leptogenesis. 

In the case I, which has the simplest form of right-handed Majorana neutrino mass matrix 
with two parameters as seen in eq. (2.28), the class S is consistent with the current neutrino 



Ml > 4.9 X 10« GeV , 



(4.8) 




(4.9) 



5 Summary 
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experimental data, if we are allowed to take a little larger value of up quark mass at the GUT 
scale. On the contrary, in the other three cases which are slightly extended to more general 
cases within 2-zero texture, having one or two new parameters as seen in eq. (2.30), (2.31) 
and (2.32), it is shown that the type 5*1 and 5*2 have the experimentally allowed regions in 
the case II, the types Ci, F3 and F4 in the case III, and the types 5*1, 5*2, Ai, Bi, Ci, F3 
and F4 in the case IV. We found that the prediction of jC/gsl depends on the types of Dirac 
neutrino and right-handed Majorana mass matrix, considerably. 

We have also calculated the branching ratios of LFV processes for the type Si and S2. 
The predicted branching ratios are well below the experimental upper bounds except r — >■ 
process for the case Si. On the other hand, for the case II, III and IV, the predicted branching 
ratios are almost same as the case I. Also, we have discussed the thermal leptogenesis in 
our model. Because in the case II, III and IV corresponding to the hierarchical case, the 
hghtest Majorana neutrino mass Mi is lighter than 4.9 x lO^GeV for all types to satisfy 
the conditions of the current neutrino experiments, it is impossible to explain the observed 
baryon asymmetry by thermal leptogenesis for the case of II, III and IV. 

In summary, we have shown that only the class S in the case I, having degenerate mass 
spectrum for the 1st and 2nd generation of right-handed Majorana neutrino, can simulta- 
neously explain the current neutrino experimental data, lepton flavor violating process and 
baryon asymmetry of the Universe. The precision measurements for neutrino mixings and 
mass-squared differences, furthermore, LFV will test if such model is realized in Nature in 
near future. 
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A Phases in neutrino mass matrices 



The complex symmetric matrices are given for the Dirac and Majorana neutrinos as follows: 

/O a 0\ /Or 0~ 



M, 



uD 



R 



r I M3 , 
1 



(A.l) 



a b c \ rrit , 
\0 c 

where each element is complex in general except for M3 and m^. The matrix Mr is trans- 
formed to the real symmetric matrix by phase matrix Pr 



R 



Mr = P*rMrP*j, 




R 





V 










I , (A.2) 
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where x — ^rg[r]. 

On the other hand, the Dirac neutrino mass matrix turns to 



M, 



vD 



where 



a 







c 


c 


\d\e 



tOLv 






















gi7. 



= X - 2 arg[a] + arg[c] , ^v^- arg[c] , 

Iv^X- arg[a] , r = arg[c] - arg[a] , 

01 = arg[a] + arg[6] - 2 arg[c] + x , 02 = - arg[a] + arg[d] + x ■ 



By using the seesaw formula, we have neutrino mass matrix 



where 



r 


c 


a| 


c 



|a||c| + |r||c||ci|e' 





+ |r||c||( 



2 

i 

r|M3 



Taking account the hierarchy of parameters \a\ ~ A^, |6| ~ A^, |c| ~ A^, \d\ ~ 
|r| ~ |a||c|/|dp with A ~ 0.2, we get 



( |aP 

lap 2|a||&|e^^^ \a\\c\ 
\ |a||c| IdlV^-^^ / l^^l^^^^ 



By using another phase matrix P', M,^ turns to 



/ 

\a\ 

\ 



2 2|a||6|e*(^i+2^2) |a||c| 
lal IcI 





z| |( 



\r\Ms ' 



where 



P' 



/e-^-^^ \ 

e*'^^ 
^ 



Therefore, the neutrino mass matrix is given as 



= P.M.P, , P, = P'*P* 







g-j(/3./+</'2) Q 

V e-^(^''-<^2) 
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Suppose the charged lepton mass matrix Mi to be real by the phase matrix Pi, 

e~ 



M, 



Ml = P^MiP; 



-lai 

























where M; is real matrix. Then, the MNS matrix Umns is given by 

Umns — OJQ*Ui, , Q — Pi,Pi , 



(A.ll) 



(A.12) 



where Ui — Pi*Oi and U^, are unitary matrices, which diagonahze M; and M^,, respectively. 
In this paper, we have parametrized 



Q 



/I \ 

e-'P 

V e-'" J 

where 

p ^ ai - - Pi + + 202 , cr ^ ai - - ^1+^^ , 

which are given in terms of the arguments of a, b, c, d and r. 

In the leptogenesis, the effective phases are 0i and 02 since we calculate in the basis of 
the real mass matrix Mr 



(A.13) 



(A.14) 



M^dMIj, = P*^M,dMIjjPr , 



(A.15) 



which is independent of phase matrix P^. The phases 0i and 02 are independent of the 
phases p and a, which appear in the MNS matrix and then, in the calculations of the lepton 
flavor violations such as — > e + 7. 



B Phases in the charged lepton mass matrix 



The complex symmetric matrix is given for the charged lepton mass matrix as follows: 

/ ai \ 

Ml = ai hi ci mi , (B.16) 
\ ci dj 

where ai, bi, ci and di are complex in general. The mass matrix turns to 





/ 





\ai\ 


,M 


( 


Mi^Mi^ Pi*MiP* = 




ai\ 


\bi\e''^' 


\ci\ mi , 






I 





\ci\ 


\di\J 


\ 











\ 


-ill J 



, (B.17) 



where there is still one phase after removing phases by the phase matrix P;. Mass eigenvalues 
and left-handed mixings are given by solving the following matrix 



m\ Ml 



( \aif \ai\\bi\e''^^ 

\ai\\bi\e-''i>^ \ai\^ + \bi\^ + |q|2 \ci\\di\ + \bi\\ci\e-"^^ \ m\ , 
\ \ai\\ci\ \ci\\di\ + \bi\\ciy^' 



+ 

|2 



(B.18) 
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Due to the hierarchy of parameters \ai\ ~ A^, ~ A^, |q| ~ A^ and \di\ ~ 1, the effect of 
the phase (pi is minor. The eigenvalue equation is approximately given as 



X 



dil'x' + {\bi\'\di\' + IqIV/I - ^IkWcilVil cos0;)a; - \ai\'\ci\'\di\' = , (B.19) 



where non-leading terms are neglected. The term including cos (pi is also a non-leading term. 
By the rephasing in eq. (B.18), the phases moves to the 1-3 and 3-1 elements as follows: 

|2 



\ai\\b, 



m\\ci\e 



i4>l 



\ai\\bi\ 
\ci\\di\ + \bi\\ci\e' 



\ci\\di\ + |&;||Q|e" 
\ci\' + W 



m 



I 1 



(B.20) 



where 



and 




(B.21) 



after neglecting non- leading terms. Therefore, the imaginary part appears only in the (1-3) 
mixing, in which the absolute value is very small compared with other mixings. In conclusion, 
the effect of the phase 0/ can be neglected in practice. 



C The effect of deviation from m2 

The following two-zero texture 



M 



/O A 0\ 

ABC 
Vo C D) 



(C.l) 



has the relations between its components and mass eigenvalues as follows: 



B + D = mi+m2 + m3, (C.2) 

BD — — A^ = mim2 + m2m3 + msTTii, (C.3) 

DA^ = -771177127713. (C.4) 

If we take B = m2(= Bq) and D — -\- mi(= Do), the remaining components can be 
obtained as [18] 



A = ./i^^^^il^(^Ao), 6-^/-^^^^^^^^-^^ + ^^^ (^Co). (C.5) 
V ms mi V ms mi 

Hereafter, we will transform mi into —mi by rephasing. Without the loss of generality, we 
can consider a small deviation e from m2 in the 2-2 component of M: 

B = m2 + e^ Bo{l + e'), (C.6) 
D = mg -mi -£ ~ £'0(1 -£), (C.7) 
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where e' — ejmii s — e/mz- Then, we obtain 



and 



= 17111712 — 11721713 + m^lTli + BD — 

~ mim2 - 17121713 + 17131111 + -Bo-Do - ^0 + BqDqs' 



Here, we can calculate 



= (C.IO) 

Gq mi 7712 mi 

In the case I (e = 0) , we needed to enlarge the range of up quark mass at the GUT scale in 
order for the class S to get the overlapped region on a-(3 plane. We, here, examine whether 
or not we can obtain the overlapped region by the effect of deviation from m2, instead of 
taking a wider range of up quark mass. 

With eqs. (C.6), (C.7), (C.8) and (C.9), the Dirac neutrino mass matrix is given as 

/ Oq \ 



M ~ 



m^, (C.ll) 



ao a^H + e') co{l + ^-^e'f'^ 

\Q co(l + ^e')'/' ^0° ) 

where we take a ~ ag and 0? ~ do because of £ -C e'. Then, the neutrino mass matrix is 
given as 

= /3 a h\ — with M« = r M3, (C.12) 

Vo h 1) ^3 Vo 1/ 

where 

h = — q; = — P = — (C 13) 

The values of a and (3 are determined by a parameter r, or equivalently h: 

a^—h, -h, (C.14) 

c c 

from which we can finally obtain 

= h ; ^-r-j T7T7 — "0 1 - 7: , (C.15) 

Co(l + Mo£m/2 "V 2mi^' ^ ^ 

/5 = — ..,0 -Po(i--—)- (C.16) 

^0 

As we can see in eqs. (C.15) and (C.16), the change of e affects both a and /3 equally, 
although we have seen that enlarging the range of up quark mass made only a decrease in 
the case I {e = 0). Therefore, we cannot arrive at the overlapped region by taking the effect 
of deviation from m2. This has been also confirmed by numerical calculation. 
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Figure 1: The black region is the experimentally allowed region predicted from a neutrino 
mass matrix with two zeros of eq. (2.1). The region predicted from the up-quark masses at 
the GUT scale for five classes S, A, B, C and F in the case I with h = 1.3. 
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Figure 2: The black region is the experimentally allowed region predicted from a neutrino 
mass matrix with two zeros of eq. (2.1). The gray region is predicted from the up-quark 
masses at the GUT scale for the type Si in the case II with h = 1.3. 
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Figure 3: The black region is the experimentally allowed region predicted from a neutrino 
mass matrix with two zeros of eq. (2.1). The gray region is predicted from the up-quark 
masses at the GUT scale for the type 5*2 in the case II with h = 1.3. 
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Figure 4: The black region is the experimentally allowed region predicted from a neutrino 
mass matrix with two zeros of eq. (2.1). The gray region is predicted from the up-quark 
masses at the GUT scale for the type Ci in the case III with h = 1.3. 
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Figure 5: The black region is the experimentally allowed region predicted from a neutrino 
mass matrix with two zeros of eq. (2.1). The gray region is predicted from the up-quark 
masses at the GUT scale for the type in the case III with h = 1.3. 
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Figure 6: The black region is the experimentally allowed region predicted from a neutrino 
mass matrix with two zeros of eq. (2.1). The gray region is predicted from the up-quark 
masses at the GUT scale for the type F4 in the case III with h = 1.3. 
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Figure 7: The black region is the experimentally allowed region predicted from a neutrino 
mass matrix with two zeros of eq. (2.1). The gray region is predicted from the up-quark 
masses at the GUT scale for the type 5*1 in the case IV with h = 1.3. 
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Figure 8: The black region is the experimentally allowed region predicted from a neutrino 
mass matrix with two zeros of eq. (2.1). The gray region is predicted from the up-quark 
masses at the GUT scale for the type 5*2 in the case IV with h = 1.3. 




Figure 9: The black region is the experimentally allowed region predicted from a neutrino 
mass matrix with two zeros of eq. (2.1). The gray region is predicted from the up-quark 
masses at the GUT scale for the type Ai in the case IV with h = 1.3. 
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Figure 10: The black region is the experimentally allowed region predicted from a neutrino 
mass matrix with two zeros of eq. (2.1). The gray region is predicted from the up-quark 
masses at the GUT scale for the type Bi in the case IV with h = 1.3. 
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Figure 11: The black region is the experimentally allowed region predicted from a neutrino 
mass matrix with two zeros of eq. (2.1). The gray region is predicted from the up-quark 
masses at the GUT scale for the type Ci in the case IV with h = 1.3. 
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Figure 12: The black region is the experimentally allowed region predicted from a neutrino 
mass matrix with two zeros of eq. (2.1). The gray region is predicted from the up-quark 
masses at the GUT scale for the type F3 in the case IV with h = 1.3. 




Figure 13: The black region is the experimentally allowed region predicted from a neutrino 
mass matrix with two zeros of eq. (2.1). The gray region is predicted from the up-quark 
masses at the GUT scale for the type F4 in the case IV with h = 1.3. 
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Figure 14: The predicted branching ratio of — > 67 process as a function of M1/2 (GeV) 
taking oq — 0, mo — 100 ~ 1000 GeV for tan /3 = 5 ~ 50 are shown. The left and right figure 
correspond to the type Si and S2, respectively. The horizontal line is experimental upper 
bound. 




Figure 15: The predicted branching ratio of r — > 67 process as a function of M1/2 (GeV) 
taking ao — 0, mo — 100 ~ 1000 GeV for tan /3 = 5 ~ 50 are shown. The left and right figure 
correspond to the type Si and S2, respectively. The horizontal line is experimental upper 
bound. 
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Figure 16: The predicted branching ratio of r — process as a function of M1/2 (GeV) 
taking gq = 0, mo = 100 ~ 1000 GeV for tan/? = 5 ~ 50 are shown. The left and right figure 
correspond to the type Si and 5*2, respectively. The horizontal line is experimental upper 
bound. 
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Figure 17: The predicted baryon asymmetry of the universe are shown as a function of 
the degeneracy AM. The region between horizontal hues is the observed value of baryon 
asymmetry ?7b [41]. The prediction is consistent with the observed value in the region of 
AM ~ 10-3. 
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